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ABSTRACT 

We analyse the chemical properties of a set of solar vicinity stars, and show that the small dispersion in abundances of 
a-elements at all ages provides evidence that the SFH has been uniform throughout the thick disk. In the context of 
long time scale infall models, we suggest that this result points either to a limited dependence of the gas accretion on the 
Galactic radius in the inner disk (R<10 kpc), or to a decoupling of the accretion history and star formation history due 
to other processes governing the ISM in the early disk, suggesting that infall cannot be a determining parameter of the 
chemical evolution at these epochs. We argue however that these results and other recent observational constraints - 
namely the lack of radial metallicity gradient and the non-evolving scale length of the thick disk - are better explained 
if the early disk is viewed as a pre-assembled gaseous system, with most of the gas settled before significant star 
formation took place - formally the equivalent of a closed-box model. In any case, these results point to a weak, or 
non-existent inside-out formation history in the thick disk, or in the first 3-5 Gyr of the formation of the Galaxy. We 
argue however that the growing importance of an external disk whose chemical properties are distinct from those of the 
inner disk would give the impression of an inside-out growth process when seen through snapshots at different epochs. 
The progressive, continuous process usually invoked may not have actually existed in the Milky Way. 
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1. Introduction 


The cosmic star formation history shows that galaxies went 
through their mos t vigorous phase of star form ation at red- 
shift z=1.5-3 (e.g Madau & Dickinson 2014|). Structures 
thought to have formed at this epoch in disk galaxies are 
classical bulges and thick disks. The stellar content of the 
nearby universe in environments with average density is 


dominated by pure d i sk galaxies dKorm endy et al. 2010 
Fisher & Drory 2010 Laurikainen et al. 2014|). Tms sug- 


gests that classical bulges represent only a minor mass com¬ 
ponent of galaxies. Thick disks, however, are found to be 
ubiquitous, an d are detected on nea rly all nearby ed ge-on 
disks galaxies (iComeron et al. 20111. Recent results (Ness 


et al. 2013a|bf Di Matteo et al. 2014a|b I have confirmed 


the view that the Milky Way is also a pure disk galaxy 


dShen et al 2010 


stantial classical 


Kunder et al. 20121, i.e., it has no sub- 
Dulge. 


Further, measurements of its star 
formation history have shown that during the peak of SF in 
the universe, the MW was a star-bursting galaxy and was 


forming its thick disk (Snaith et al 


bursting 

r]|2014ait 


b Lehnert et al. 


20141. The structur al parameters of this population, with 
a short scale leng th (Bensby et al. |2011[[Cheng et al.|2012a| 


Bovy et al.|2012 1 suggest it constitutes a major part of the 


mass in the Milky Way’s central region (Ness et al. 2013a,b, 
Di Matteo et al. 2014a,b, Haywood et al. in preparation). 


Send offprint requests to: M. Haywood, e-mail: 

Misha.HaywoodSobspm.fr 


This emerging scheme suggests that the formation of 
the thick disk is crucial phase in the formation of disks (see 


Lehnert et al. 20141, both because this population could 
represent a substantial fraction of the mass of galaxies and 
because it may account for the main increase in the metal 
enrichment of galaxies since the formation of the first stars. 
These views are opposite to the standard chemical evolu¬ 


tion scenario based on long term infall (e.g Chiappini et 


al. 1997 Fraternali 20131 and the recognition that the 


solar vicinity contains too few low metallicity dwarf stars 
compared to the “closed-box” model. In other words, ob¬ 
servations of the solar vicinity would suggest that the en¬ 
richment rate has been faster than the star formation rate. 
Infall models have been introduced in order to solve this 
problem by limiting the dilution of ejected metals to a disk 
parsimoniously provisioned with gas, and by boosting the 
increase of the metallicity at early times without creating 
too many metal-poor stars. Hence, these scenarios leave no 
room for a substantial component of intermediate metallic¬ 
ity in the Milky Way. This credo has been questioned by a 


new scheme pres ented in H aywood e t al. (20131, Snaith et 
( 2014a|b I and Haywood J 2014a|bp l 


Another key feature of the standard chemical evolution 
scenario is to make the accretion of gas radially dependent, 
with slower accretion outwards, in order to mimic an inside- 
out formation. Inner regions form their stars first, because 
gas density increases faster, while outer regions grow at 
later times and at a slower pace. Because of the coupling 
assumed between accretion and SFR in these models, a nat- 
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2. Data 
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We use the sample of|Ad i bekya n et al. (20121 already anal¬ 
ysed in Haywood et ^ ( |2013p . I'he sample contains 363 
stars with age determinations, t o which we added further 
the 36 stars from the sample of Nissen & Schuster] (20101 
for which we could determine ages. Ages for this last sam- 
ple were determined in the same manner as those from 
Adibekyan et al. (20121 sample in Haywood et al. (2013). 
However, as reported by Schuster et al. ~ f2012p , the etfec- 
tive temperature of the stars in Nissen &: Schuster may have 
been underestimated. We checked this statement by com¬ 
paring the 11 Hipparcos stars in common between the two 


samples. All these objects have a temperature lower in Nis- 


Fig. 1. The [a/Fe]-[Fe/H] distribution at the solar vicinity for 
our sample. The colors (indicated in the legend) and size (larger 
implies older) of the symbols Indicate the age of the stars. The 
stars below and to the left of the line are either accreted stars 


sen & Schuster (2010) than in Adibekyan et al. (2012), with 
a mean otfse t of -75K, confirming the statement ot |Schuster| 
et al.| (2012|). Following their recommendation, we added 
5 k to the ~ 


+ 100k. to the Nissen & Schuster (2010) values. 


according to Nissen & Schuster (20101 (for those objects which 
have [Fe/H]<-0.7 dex), or metal-poor thin disk stars. They are 
represented by empty circles in Fig. 2, and removed from the 
age-alpha distribution in Fig. 3. 


In Haywood et al. (2013) the spectroscopic scale of 
Adibekyan et al. (|2012| was compared to photometric tem¬ 
peratures using the Teff-V-K calibration of [Casagrande et] 
al. (2010), from which we deduced that the two scales are 
on by 56K. The impact of this offset on stellar ages of 
such differen ce in temperatures was studied in |Haywoo'd| 
et al. (2013) and found to be limited (see their Fig. 

Ages were determined usin g the bayesian method of |j0r-| 
gensen & Lindegren (2005), and implies the derivation of 
a probability distribution function (or PDF) for each star. 
We adopted the Yonsei-Yale set of isochrones, with atmo- 
sph eric p arameters taken from the study of [Adibekyan et| 
al. (2012), including mean [a/Fe] values. The uncertainties 
on the ages are derived fro m the PDF as described in |j0r-| 
gensen & Lindegren (2005), and the values (see Fig. 3) are 
for a confidence level of 68%, corresponding to the ±1 cr 
Gaussian errors. Random uncertainties on age reflect the 
uncertaintie s on the atmospheric parameters. In [Haywoodj 
et al. (2013), we assumed 50K uncertainty on effective tem- 


Fig. 2. The age-[a/Fe] distribution of stars in Fig. 1. Metal- 
poor thin disk stars and accreted stars, as selected from Fig. 1 
(below and to the left of the limits in gray), are shown as white 
dots. 


ural signature of such scheme would therefore be that star 
formation histories vary with distance from the galactic cen¬ 
ter. Direct measurement of the star formation history out¬ 
side the solar vicinity is not feasible currently, however we 
know that a given SFH leaves specific signatures on the ra¬ 
dial and temporal dependence of the chemical abund ances 
in stars ( Haywood et al.|[2013 jSnaith et al.||2014a|b| ) Since 
thick disk stars observed in the solar vicinity were born 
over a wide range of galactocentric distances (see section 
[^, their chemical properties should therefore reflect the 
variety in the radial star formation histories of the thick 
disk, if any. 


In the following, we show that these chemical properties 
are in fact remarkably homogeneous at all ages in the thick 
disk, leaving little room for some significant variation in 
the star formation history of this population as a function 
of radius. The outline of the paper is as follows: in the 
following section, we briefly describe our data, in Section 
we describe the model and present our results, and discuss 
them in Section [H 


peratures, noting that the formal uncertainties quoted by 
Adibekyan et al. (2012) are sometimes much smaller than 
this, 0.1 dex on metallicities, and 0.1 magnitude on abso¬ 
lute magnit udes. The total uncer tainties in the ages, as 
described in Haywood et al. (2013), are of-order 1 Gyr for 
the random error (with a mean uncertainty of 1.1 Gyr), 
and a systematic error due to stellar physics of 1 to 1.5 
Gyr. The systematic uncertainties were estimated by com¬ 
paring the ages determined using different sets of isochrones 
(see Haywood et al. 2013 for details) Given the dispersion 
of points on Fig. 3, of the order of 1 Gyr, we deemed our 
error estimates on age as realistic. 

The [Q:/Fe]-[Fe/H] distribution of the whole sample is 
shown on Fig. 1. 

The stars below the separating line are considered either 
as metal- poor thin disk ob jects (see discussion for th ese 
objects in Haywood (2008) and Haywood et al. ([2013 ) for 
the status of these objects and Snaith et al. p014a|b) for 
a description of their chemi cal evolution), or as accr eted 
stars, following the study by Nissen & Schuster (2010). 

In Fig. 2, which shows the age-[a/Fe] distribution (the 
color now coding the metallicity of each star), these objects 
are represented by empty circles. The final distribution of 
Age-[a/Fe], after removing accreted and metal-poor thin 
disk stars, is plotted in Fig. The sequences of the thin 
and the thick disks appear clearly on the plot, with a sharp 
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Haywood et al.: Clues to the formation of the Milky Way’s thick disk 



Fig. 3. The age-[a/Fe] distribution at the solar vicinity after 
having removed metal-poor thin disk stars and accreted stars, 
as selected from Fig. 1. Note the remarkable small dispersion in 
both the thin and thick disk sequences of stellar age and [a/Fe]. 
Errors on age determinations are described in the text. On 
[a/Fe], 0.02 de x error is taken as repr esentative of the internal 
errors given in 


Adibekyan et al. (20121. 


change o f slope between the tw o at 
cussed in 


Haywood et al. (20131. 


.^8 Gyr, as first dis- 


No w more clearly demonstrated than in [Haywood et al.j 
(2013||, a number of objects with a-enhanced abundance 
ratios have age lower than the sequence of thick disk ob¬ 
jects (age<8 -9 Gyr at [Q!/Fe[^0.2-0 .3 dex). As we already 
discussed in Haywood et al. (2013), these objects have no 
peculiarities.^ We did not note at that tim e however the 
dedicated study of Fuhrmann et al. (2012) on HIP 38908 


(HR 3138). jFuhrmann et al. (|2012p found similar “young” 
age to what we d erive tor this object: we derive an age of 


5At.2^ Gyr, while 


Fuhrmann et al. 


(|2012| give Gyr. 

that radial velocity ex- 


They commented further on the tact 
eludes the possibility that HR 3138 is a blue straggler with 
a white dwarf. The authors could not find any specific rea¬ 
son for the young age, but suggested that either the star 
had accreted gas from interstellar material or that the en¬ 
hanced effective temperature and luminosity resulted from 
a merger with a companioi|^ 


In Fig. 4, we plot [a/Fe] as a function of the pericenter 
distance of stars. The purpose of this plot is to illustrate 
that the pericenters of thick disk stars span the whole range 
of radii from 8 kpc to the innermost regions. Hence, un¬ 
less the orbital characteristics of these stars were seriously 
affected by dynamical evolution, it can be assumed that 
they represent a reasonable sampling of the thick disk over 
a wide, perhaps even the whole, range of radial distances 
from the Galactic center. 


^ The other stars in our sample with similar characteristics - 
HIP 69780, 77637, 51191, 54469, 77610, 109822- have no known 
peculiarity. Most (HIP 69780, HIP 77637, HIP 51191, HIP 54469 
and HIP 109822) have clear asymmetric drift and overall hot 
kinematics, confirming that they should belong to the thick disk 
or halo, in spite of their age. These objects are not considered 
further in the present study. 



d 


0 2 4 6 8 10 

Rmin [kpc] 


Fig. 4. Mean alpha abundance as a function of pericenter 
distance. Color (indicated in the legend) and size of the symbols 
indicate the age (larger symbols implies older ages) of the stars. 
Thick disk stars are essentially those stars with [a/Fe]>0.1 dex. 
Grey points are for the whole sample. 


3. Results: Modelling the [a/FeJ-age relation 

3.1. The model and chemical tracks 


The model is a closed- box and is discussed at length in 
Snaith et al. (2014a|bl. Note that the effects described 
here are expected to be of similar, or even of higher ampli¬ 
tudes with standard infall models. The reason being that 
since the ISM is less important in these models, they are 
more sensitive to the level of metal injection in the ISM 
than closed-box models, which have greater inertia in this 
respect. A clear illustration of the effect of the variation of 
the SFHs on increasing t he dispersion of alpha a bundances 
at a given age is given in Minchev et al. (2013), their Fig. 
7. 

In the following, we use the n omenclature of the disks 
as given in [Haywood et al. (20131. That is, we consider the 
thick and thin disk up to radial distances of 10 kpc to form 
the inner disk, while stars beyond this limit are members 
of the outer disk. We sketched a scenario for the formation 


of these two different structures in Haywood et al. (2013) 
and the chemical evolutions of both are described in lSnaithl 
eTHl (120140. 


In short, the closed-box model assumes that: (1) all 
the gas is in the system at the beginning. We take the 
closed-box model as the approximation of the case where 
most accretion in the inner disc (R<10 kpc) occurred early. 
(2) The ISM is initially devoid of metals and remains well 
mixed during its whole evolution. (3) No coupling between 
th e star format i on rate a nd the gas content is assumed. 
In Snaith et al. |( |2014a|b| ), we determined the SFH by fit- 
ting our model to the Age-[Si/Fe] data. Here we assume 
different star formation histories from which we derive the 
Age-[Si/Fe] tracks which are compared to the data. 

Th e ingredients of the model are: standard theoretical 
yields (Iwamoto et al. 1 999 Nomoto et al. 2006 Karakas 
2010 1, an IMF (Kroupa 200ip independent of time, and a 
stellar mass-lifetime relation dependent on the metallicity 


(Raiteri et al. 1996). No instanta neous recycling appr oxi¬ 
mation has been implemented, see Snaith et al. (2014b I for 
details. 
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3.2. Comparison with the data 


Models with radially dependent accretion time scales have 
a range of star formation histories that vary from having a 
peak at early times in the inner disk, to constant, or even 


increasing 

SFH in the solar vicinity or beyond (e.g. Naab 

& Ostriker 

[2006 Minchev et al.|2013 Kubryk et ai.[|2014|). 

Hence, in . 
30 M 0 .pc“ 
10 M 0 .pc“ 

VLinchev et al. the Sh'R at R=4kpc peaks at ~ 
^.Gyr“^ in the first Gyr, decreasing to less than 
■^.Gyr“^ at the present times, while it is only 


slightly decreasing (by a factor of ^ 2 ) at the solar radius, 
and very nearly constant at 12kpc. In Naab & Ostriker, the 
SFH at 4 kpc peaks a later time (6 Gyr), reaching similar 
values (20-30 M 0 .pc“^.Gyr“^), and decreasing to about 10 
M 0 .pc“^.Gyr“^ at the present epoch. The peak of the 
SFR at the solar vicinity is reached only at ~11 Gyr and 
thus is mostly increasing throughout the evolution of the 
MW. Independently of any particular SFH, these models 
show that the SFR at early times varies by large amounts 
depending on the radius: there is a factor of about 5 to 10 
at 2-4 Gyr between the SFR in the inner regions (2-4 kpc) 
and the outer regions (at about 8-10 kpc from the Galactic 
center). 

Because the radial variation of the SFHs can vary widely 
between one model to the other, we do not rely on any par¬ 
ticular model. We have chosen a set of SFHs which we take 
as representative of the general phenomenology of radially 
dependent SFHs in disks, and that are shown in Fig. 
They are taken as illustrating the range of variation that 
one can expect from adopting an infall time-scale and re¬ 
lated SFH varying with galactic radius: increasing SFH, 
constant SFH, decreasing SFH, and reflect the 5-10 ampli¬ 
tude of the SFR in the first few Gyr as noted above. 


Fig. | 6 ] shows the corresponding chemical tracks gener¬ 
ated from these SFHs over-plotted on the data. As can be 
seen, the range of alpha abundances spanned by the mod¬ 
els at any given age is large, bracketing the observed [Si/Fe] 
values. As expected, the SFH that fits best the observed 
age-[Si/Fe] distribution is one with a significant peak of star 
formation at e arly times, similar to the one deduced by fit¬ 
ting the data in Snaith et al. (2014a I. If the peak in the SFH 
occurs too early (< 2 Gyr), the model generates too much 
iron early on and the chemical track is too low compared 
to the data. However, if the SFH is flat (constant SFR), 
the knee at ^ 8 Gyr is not reproduced, with the chemical 
track being off the data by about 0.1 dex over this region. 
In the case where the SFR starts low and increases to reach 
a maximum in the last few Gyrs as suggested by the model 
of [Naab & Ostriker (2006|for the solar vicinity, the model 
is far from the data. 


As shown in the previous section, the pericenter dis¬ 
tances of thick disk stars span a large radial range from less 
than 2 kpc from the Galactic center up to the solar circle 
(Fig-i- Hence, we do expect thick disk stars to sample the 
whole thick disk (R< 10 kpc) and their alpha abundance 
ratio to reflect star formation histories at a variety of radii. 
Obviously, the expected dispersion is not observed. The 
thick and thin disk sequences are very well defined and at 
all ages, the dispersion in abundance ratio is small in both 
phases. In fact, the observed dispersion at a given [Si/Fe] 
is fully compatibl e with the uncertainty of about 1 Gyr on 
the derived ages ( Haywood et al.|[2013 l, suggesting the in¬ 
trinsic dispersion in abundance is essentially negligible. In 
particular, we would expect the transition from the thick 
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Fig. 5. The set of star formation histories used in this analysis. 


to the thin disk at ~8 Gyr to be much less well defined 
than it is. Had the thick disk SFHs been different in the 
inner and outer regions, the alpha abundance ratios at the 
end of the thick disk phase would show a level of dispersion 
in their a bundance ratios which is not observed. As com¬ 
mented in Haywood et al. (2013), the data strongly suggest 
that the thick disk phase ended by leaving a disk of uniform 
properties from which the thin disk started to form. 

While the previous discussion was meant to be illustra¬ 
tive and useful for comparing with SFH commonly used in 
models, the most pertinent question is really: How much 
variation in the SFH is really allowed by the data? This 
can be estimated by using the bootstrap SFHs generated 


Snaith et al. (2014b I to evaluate the uncertainty in the 


mean Sh'H that tits the data. The bootstrap SFHs give 
us some kind of “cosmic variation” that is allowed by the 
constraints provided by the data. What we seek here is 
not the uncertainty in the SFR at a given time, but how 
much the SFH can be systematically below or above the 
best fit SFH over a given time interval. In order to do so, 
we summarize the information given by each of the boot¬ 
strapped SFHs by measuring the ratio of the stellar mass 
formed during the thick disk phase (13-8 Gyr) to the to¬ 
tal stellar mass formed (13-0 Gyr), for each solution. The 
mean and dispersion of these percentages are 56 and 2.5%, 
with values that vary, at most, between 50 and 61%. Note 
that the derived uncertainty in the SFR at a given age is 
much l arger than 2.5% (see t he error bars in the derived 
SFH in [Snaith et al. (2014a|b)), but the overall variation of 
the SFR, as measured by the percentage of mass created in 
the age range 13-8 Gyr relative to the total stellar mass, is 
much better constrained. Gomparatively, the stellar mass 
produced between ages 13 and 8 Gyr - during the thick 
disk phase - relative to the total stellar mass in each of the 
model SFHs (Fig. 5) is 11, 26, 33 and 63%. Hence, any 
significant varia tions from the mean SFH derived in [Snaith] 
et al. ( [2014a|b[ | are ruled out. 


4. Discussion 

Summarizing the previous section, the observed dispersion 
in the age-[Si/Fe[ distribution essentially reflects the uncer¬ 
tainties in the derived ages, although the thick disk stars 
originate from all over the disk. This implies that the chem¬ 
ical evolution in the thick disk phase must have been an 
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Age [Gyr] 


Fig. 6. The age vs [Si/Fe] chemical evolution tracks corre¬ 
sponding to the SFHs given in Fig. 5, compared to the stellar 
age vs [Si/Fe] data. The colors of each line correspond to those 
in Fig. 5. Thick part of the curves correspond to the thick disk 
phase. 


extremely homogeneous process, incompatible with signifi¬ 
cant spatial variations of the SFR at these epochs. These is 
in clear contradiction with a number of models, which ex- 
plicitly predict a very significant dispersion, see for example 


Minchev et al. ( 2013[ ), or Kubryk et al. (20141, which both 


have a range of [a/Fe] larger than 0.2 dex at the end of the 
thick disk phase. Note however that this is likely to be true 
for all models which, motivated by the inside-out scenario, 
build there disk through a radially dependent SFH. 


Another, perhaps more natural explanation, could ex- 
plain t he seemin gly u niform SFH in the thick disk. In jHay-j 
wood ( 2014a|b ) and Snaith et al. ( 2014a|b |, we provided 


a number of arguments to propose a new scenario for the 
chemical evolution of the Milky Way disk. We argued that 
the structural parameters and the measurement of the SFH 
of the Milky Way provide new evidence that the thick disk is 
massive, and that this is not expected in standard chemical 
scenarios (see the Introduction). The only possible way to 
generate sufficient numbers of stars at intermediate metal- 
licities (which, we argued, are mostly confined to the inner 
regions, as o ne would expect for a population wi t h 2 kp c 
scale length, |Bensby et al.| (|2011|), |Cheng et al.| (|2012aD, 
Bovy et al. ( 2012p ), is that the gaseous disk was already 


massive at early times. From the point of view of chemical 
evolution, the model that is the most likely to represent 
this situation is a closed-box model. In such case, because 
of the short scale l ength of the disk, th e relative large scale 
of the turbulence (|Lehnert et al.||2014|), and the absence of 
radially dependent accretion time scale, less differentiation 
in the SFHs is expected to take place in the inner and outer 
regions of the thick disk. 


Nidever et al. (20141 confirm this conclusion on the basis 
of the analysis of the APOGEE data and a chemical evo¬ 
lution model. Their model includes infall with a long time 
scale, but no radial dependence, hence, no induced varia¬ 
tion in the SFH due to some kind of inside-out formation. 
They assessed however possible variations in the SFH by 
assuming different star formation efficiency, and deduced 
that some homogeneity in the SFH can be deduced from 
the limited dispersion observed in the [a/Fe]-[Fe/H] distri¬ 
bution in t he thick disk regime s amp led by the APOGEE 
data. As in|Haywood et al. (2013) and Lehnert et al.|(2014), 
they note that this implies that the thick disk has formed 
from a turbulent, well-mixed gaseous disk. 


4.1. The radially uniform SFH of the thick disk 

There are a number of important consequences of this re¬ 
sult. In chemical evolution models based on long-time scale 


al. 


2013 Kubryk et ml 2014p , the SFH depards on the 


accretion of gas (e.g. jChiappini et ah |[T997l jMinchev et 


with caution however, because a s shown recently (Hopkins 
et al.l 2014 Woods et al. 20141, the SFH does not nec- 


iniall history through some form of relation between the 
star-formation intensity and the gas-surface density (i.e. a 
Schmidt-Kennicutt law). The infall law is chosen in order to 
reproduce the phenomenology of inside-out galaxy forma¬ 
tion models, and one of their possible signatures, chemical 
radial gradients. Our results contradict a fundamental as¬ 
pect of these scenarios: since the SFH is measured to be the 
same across the (thick) disk, it may likely imply that there 
is no radial dependence of the accretion of gas onto the in¬ 
ner disk (R<10 kpc). Although such simplistic inferences 
are similar to those assuming direct connection between in¬ 
fall and star formation intensity, they should be interpreted 


4.2. On the inside-out galaxy formation model 

Other evidence supports a uniform SFH for the thick disk. 
First, no radi al metallicity gradie nts have been found in 
the thick disk ( Cheng et al.|201^ I. This is expected if the 
SFH is radially undifferentiated. Note that the tight cor- 
relation found between age and chemical parameters (see 
and previous section) in the thick 
ack of gradient is not a consequence 


Haywood et al. (20131 
disk implies that the 


of mixing in the stellar disk, but of a well mixed I SM at a 
given epoch in the thick disk (Haywood et al.|2013 1. Chem¬ 
ical properties of the thick disk e volved homogeneously over 
time. We already com mented in Haywood et al. (20131 and 
Lehnert et al.|(2014l that this was most probably a conse- 


quence of the high turbulence in the gaseous disk at these 
epochs. Turbulence induces sufficient mixing to ex plain the 
high homogeneity of chemical species in the ISM (Feng & 


Krumholz 2014). Stars formed out of this well-mixed gas 


essarily correlate with the accretion history. Feedback is 
expected to play an important role by reducing star for¬ 
mation rate and efficiency thereby producing a reservoir of 
gas which can affect the long-term fueling of the SFH of a 
galaxy. Hence, infall could be radially dependent, but not 
directly correlated to the SFR, implying that it cannot be 
a determining parameter of the chemical evolution in the 
inner disk. 


in a geometrically thick layer. 

Second, no dependence of the radial scale length of the 
thick disk as a function of the g -abundance ratio s has b een 
detected (see our comments in Haywood et al. ( 2013|) on 
the Fig. 5 of Bovy et al. (20121). Haywood et al.j (|i013| 
demonstrated that alpha abundances correlate well with 
age in the thick disk phase (see also the previous section). 
Hence, this implies that no significant dependence of the 
scale length with age, or an inside-out age trend during the 
formation of the thick disk, is detected. This is in accor- 
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dance with the fact that no significant variation of the SFH 
is imprinted in the age-a relation, which we would expect 
if the youngest populations had the longest scale lengths. 
This is al so in good accord with the general scheme pre¬ 
sented in Haywood et al. (20131 where we discussed the 
inside-out galaxy formation scenario. 

Note that the Haywood et al. study has been cited as 
claiming evidence that the MW did not form inside-out. 
We caution that this is an oversimplification of the view 
presented in that paper. What we claimed is that there 
is no such structure as a single disk formed smoothly and 
continuously from the inside-out and generating smooth ra¬ 
dial chemical gradi ents. Indeed, there are two disks in our 
nomenclature (see Haywood et al. (2013)): the inner one 
(R<9-10kpc), composed of the thick and th e thin, with 
differe nt scale length^~2 and 3.6 kpc, see Bovy et al. 
(|2012||), but with no evidence of an increase in the scale 
length as each of these disks grew. And there is the outer 
disk, which started to form before the inner thin disk, 10 
Gyr ago, and which has a long scale length (^4 kpc or 
more). Hence, the general trend at early times (> 8 Gyr) 
would be mostly imprinted by the formation of the mas¬ 
sive, and concentrated, thick disk, followed by the outer 
(but less dense) disk which has a much longer scale length. 
Hence, snapshots of such a galaxy at different epochs would 
certainly give the impression of an inside-out growth of its 
disk, although the overall sequence is more complex and is 
not a smooth process. What we did claim in Haywood et 
al. and now argue further is that this overall scenario is in 
contradiction with the simplest picture of inside-out galaxy 
formation as we now discuss. 


4.3. The formation of disks 

Within the context of these results, models of formation of 
MW-like galaxies may need some, perhaps significant, re¬ 
visions. Models where the radial dependent gas accretion 
drives a sequence of star formation which is inside-out dur¬ 
ing the formation of the thick disk are now likely ruled out. 
In addition, models, at least during the phase of thick disk 
growth, need to be such t hat the exponential disk s cale 
length does not increase. 


Ferguson & Glarke 


(2001) ex- 


^ Note that 
[Fe/H]>-0.1 
for the thin 
abundance, 
to belong to 
shorter than 
in Haywood 


the value of 3.6 kpc is a represent ative mean, but at 
dex the scale length measured by Bovy et al.|( 20121 
disk is a complex function of metallicity and alpha 
For instance, metal-rich stars, which are more likely 
the inner thin disk have a scale length significantly 
the value of 3.6 kpc. This complex issue is discussed 
et al. (2015, in preparation). 


existing disk rich in gas. Interestingly, Ferguson & Glarke 
predict that, in their pre-assembled disk, “half of the stellar 
mass is in place more than 6 Gyr ago”, in good ag reement 
with the SFH measured by Snaith et al. ( 2014a|b I. 

The advantage of viscous dissipation is that it natu rally 
forms exponential dis ks with a large span in radius (Fer¬ 
guson & Glarke 2001) and because the viscous time scale 
is similar to that of tne star formation timescak^ it helps 
to alleviate the problem of forming massive bulge s which 
plague simulations witho ut strong stellar feedback ([Dutton 


& van den Bosch 2009). Viscosity both moves gas out- 
wards, extending the disk, and prevents gas from accumu¬ 
lating into the center of the potential. Since the disk scale 
length does not increase when both viscosity is important, 
and the accretion and star formation timescales are similar, 
gas is distributed with the same radial dependence as the 
previous generation of stars. 

Thus viscosity is able to prevent the thick disk from de¬ 
veloping a radial metallicity gradient or increasing its scale 
length with time or with increasing metallicity. 

(1^1 

h 


We note that the analysis of Vollmer & Leroy 
suggests that star formation in local disk galaxies wit 


. stel- 


plored how the scale length of disks is affected during mfall 
when the time scales for star formation and viscous gas 
redistribution are similar (within a factor of a few). They 
show that, if the main infall phase precedes the onset of star 
formation (the equivalent of a closed box model), the scale 
length is invariant with time, because the scale length is 
determined by the mean specific angular momentum of the 
gas already settled and viscosity redistributes the angular 
momentum such that the scale length is invariant even as 
the total extent of the disk grows. On the contrary, if the in¬ 
fall time scale is similar to the lifetime of the disk, the scale 
length evolves according to the specific angular momentum 
of the accreted gas. Hence, at first order, a non-evolving 
scale length in the thick disk is a natural outcome of a pre- 


lar mass M* ^10^° Mq may also be sustained by the viscous 
dissipation of the (outer) HI extended disks. Viscous dis¬ 
sipation allows the gas to move inwards and distributes it 
radially within the disk at rates sufficient to support the 
observed rate and distribution of star formation. Their 
analysis of local disk galaxies and our results imply that 
viscous dissipation and gas redistribution may be generi- 
cally important to the growth of disks. 

Gas accretion, regulated by the intensity of on-going 
star formation, may help to alleviate both the problem of 
needing the gas accretion to follow the radial distribution of 
the star formation an d of maintaining a cons tant disk scale 
length . For example, Marinacci et al. (2011) and Marasco 
et al. (2012) have suggested gas driven into the hot coro- 
nae by recent star formation in the disk may cause the halo 
gas to cool. This cool gas can then accrete onto the galaxy. 
Since the rate at which mass and momentum are injected in 
the halo gas controls how much gas is ultimately accreted, 
and these injection rates depend on the star formation in¬ 
tensity of the disk, generates a self-regulating cycle whereby 
more intense star formation leads to more intense accretion, 
resulting in further star formation, and so on. 

In such a model, unlike the recent models of gas a ccre- 
tion from the cosmic web (e.g. Danovich et al. 2014), the 
gas accretion rate would follow the radial intensity of on¬ 
going disk star formation and, perhaps in conjunction with 
viscosity, maintain a constant disk scale length. 


The results of the study by van Dokkum et al. (2013 (fits 
well in this scheme, first because of a similar mass growth 


^ The viscosity is given hy u = vturb kurb, where vturb is the 
turbulent velocity of the largest energy injection scale, and Uurb 
is this largest scale. The viscous time scale is simply, ti^ = R^/i' 
(Vollmer & Leroy 20111. Scaling this relation to observations 
of z=2-3 galaxies and values consistent with the thick disk of 
the M W ( Bovy et al.||201^ Haywood et al.|2013 Lehnert et al. 
20141, suggests, u = 0.05 (vturb/50 km s“^) (hurb/lOOO pc) pc^ 
yr“ and t^ = 6x10^ (R/1.8 kpc)^ yr. The global rate at which 
stars form is approximately one dynamical time of the thick disk, 
tdyn= 6x10^ (R/1.8 kpc) (vrot/200A:ms“^)^ yr. These crude es¬ 
timates suggest that plausibly, t„~tdyn or similar to the star 
formation time scale of the disk, and thus viscosity may be im¬ 
portant for redistributing the gas. 
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found on progenitors of the M ilky Way (see Snaith et al. 
(2014aI, Lehnert et al. (2014|), but also by showing that 
their galaxies do not increase their size through inside-out 
growth between z=2.5 and 1 , but through self-simil ar mass 
growt h at all radii ( see al so Morishita et al. (20151). Sim¬ 
ilarly, Nelson et al. ( 2012[ | find that ga laxies begin to form 
inside-out after z=l (see comment by van Dokkum et al. 
T^). This is precisely what is predicted by viscous re- 
istribution of the gas (Ferguson & Clarke 20011, and in 
agreement with our findings here. 


5. Conclusions 


We have studied the impact of varying the SFH on the 
age-[a/Fe] relation of disk stars. We show that a radially 
dependent growth of the disk would give rise, through a 
radial variation of the SFHs, to a dispersion in this relation 
which is not observed. Star formation in the thick disk was 
a general, large scale process, that led to a homogeneous 
chemical evolution throughout the inner disk (R<10kpc). 
Consequences of these results for chemical evolution mod¬ 
els have been explored. If infall intervenes on a long time 
scale, as advocated in standard chemical evolution model, 
it implies either that infall is not varying with Galactocen- 
tric distance in the inner disk, or that infall is decoupled 
from the SFH. In any case, it implies that infall is not a 
determinin g parameter of the c h emical evolution of th e in¬ 
ner disk. In Haywood|(|2014a|b|),|Snaith et al.|(|2014a|b|), we 
have shown that the characteristics of the thick disk favor a 
scenario where stars started to form from a disk rich in gas 
(most accretion occurred before significant star formation 
took place), which we model as a closed-box. We argued 
that the present results find better justification within this 
new context. It is reasonable to expect that a uniform SFH 
is more likely to occur in a disk with a short scale length and 
given the high level of turbulence in disks at high redshift, 
and with no radially dependent infall time scale. Moreover, 
differentiated SFHs are expected to produce both radial 
chemical gradients and an increase in the scale lengths of 
the disk, expected in inside-out formation scenarios, while 
none are s een in the observations. The study by |Ferguson| 
& Clarke (20011 is qualitatively in accord with our results: 
they argue that, under the assumption of similar time scale 
for star formation and viscosity, in a pre-assembled disk, 
no significant variation of the scale length is expected, and 
therefore no inside-out growth of the disk. 
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